Casimir Interactions for Anisotropic Magnetodielectric Metamaterials 
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We extend our previous work [Phys. Rev. Lett. 100, 183602 (2008)] on the generalization of the 
Casimir-Lifshitz theory to treat anisotropic magnetodielectric media, focusing on the forces between 
metals and magnetodielectric metamaterials and on the possibility of inferring magnetic effects by 
measurements of these forces. We present results for metamaterials including structures with uni- 
axial and biaxial magnetodielectric anisotropics, as well as for structures with isolated metallic or 
dielectric properties that we describe in terms of filling factors and a Maxwell Garnett approxima- 
tion. The elimination or reduction of Casimir "stiction" by appropriate engineering of metallic-based 
metamaterials, or the indirect detection of magnetic contributions, appear from the examples con- 
sidered to be very challenging, as small background Drude contributions to the permittivity act to 
enhance attraction over repulsion, as does magnetic dissipation. In dielectric-based metamaterials 
the magnetic properties of polaritonic crystals, for instance, appear to be too weak for repulsion 
to overcome attraction. We also discuss Casimir-Polder experiments, that might provide another 
possibility for the detection of magnetic effects. 

PACS numbers: 12.20.-m, 78.20.Ci, 81.05.Zx 



I. INTRODUCTION 

The last decade has witnessed a hiige activity in the 
development of metamaterials (MMs) [l| , boosted by the 
possibihty that such engineered media may give rise to 
novel optical properties at selected frequency ranges, in- 
cluding negative refraction [3], perfect lensing [3|, and 
cloaking [j, among others. Such striking phenomena, 
inaccessible with natural materials, are all possible due 
to the significant magnetic activity built into metamate- 
rials, starting at microwave frequencies and going all the 
way up to the optical range. Generally speaking, meta- 
materials arc made of micro- and nanostructures care- 
fully designed to collectively endow them with a partic- 
ular electromagnetic property. It is generally desirable 
that these structures should be smaller than the wave- 
length of the incident radiation, so that they are seen by 
the incoming waves as artificial "atoms". This fact of- 
ten allows the use of an effective medium approximation 
to describe metamaterials in terms of an effective elec- 
tric permittivitiy tensor e{Lo) and an effective magnetic 
permeability tensor h{uj), which incorporate the typical 
optical anisotropy of metamaterials. 

Recent years have also witnessed an increased interest 
in Casimir physics [a, la] thanks to improved precision 
measurements [3] of the force between material objects 
separated by micron and sub- micron gaps. Quantum vac- 
uum fluctuations are modifled by the presence of mate- 
rial boundaries, and this typically results in an attrac- 
tive Casimir force that depends sensitively on the shape 
and the optical properties of the boundaries. While the 
Casimir force offers new possibilities for nanotechnology, 
such as actuation mediated by the quantum vacuum, it 
also presents some challenges, as micro- and nanoelec- 
tromechanical systems (MEMS and NEMS) may stick 
together and cease to work due to the attractive na- 
ture of van der Waals and Casimir forces. A strongly 



suppressed Casimir attraction, or even repulsive Casimir 
forces, would provide an anti- "stiction" eff'ect. Repulsive 
Casimir forces between two objects 1 and 2, immersed 
in a background medium 3, may come in a variety of 
ways. One possibility involves non-magnetic media only, 
for which repulsion happens when the electric permittiv- 
ities evaluated at imaginary frequencies satisfy the re- 
lation €i{iuj) < e3{iuj) < e2(iu!) [8|. Another possibil- 
ity, first predicted by Boyer [9|, involves magnetodielec- 
tric media: there is a repulsive force when a perfectly 
conducting plate is placed near a perfectly permeable 
one with vacuum in between. Some years later it was 
shown that Casimir repulsion can also occur between 
real {i.e., non-ideal) magnetodielectric media, as long one 
medium is mainly electric and the other one is mainly 
magnetic [10[. However, this possibility has been con- 
sidered unphysical [111 ] , as naturally occurring materials, 
even strong magnets at low frequencies [i2l , do not show 
significant magnetic response at near-infrared and opti- 
cal frequencies, which has been assumed as a prerequisite 
for repulsion between Casimir plates separated by typi- 
cal experimentally relevant distances of d = 0.1 — 1/im. 
On the other hand, recent developments in nanofabri- 
cation have resulted in metamaterials with magnetic re- 
sponse in the visible range of the electromagnetic spec- 
trum 1^14^ 15j. fueling the hope for Casimir repulsion 



[la . \Vn . [Tsl . IT9| . The expectation is that, by tuning this 
magnetic response to the right frequency range and mak- 
ing it strong enough, one could produce an experimen- 
tally measurable Casimir repulsion between, say, a MM 
slab and a thin metallic plate, or at least a significantly 
reduced attraction. 

Unfortunately, this is easier said than done. The major 
issue is that the Casimir force between real dispersive ma- 
terials is a broadband frequency phenomenon, as shown 
by the Lifshitz formula expressing the force between two 
semispaces as an integral over all (imaginary) frequencies 



with an exponential cut-ofF c/d [20] ■ For Casimir repul- 
sion purposes, this requires a magnetic response strong 
enough to dominate the electric response of the material 
in a broad range of frequencies, which typically is not 
the case for the magnetic resonances present in metama- 
terials. In addition, several metamaterials have metallic 
inclusions that produce a low-frequency Drude-like be- 
havior in e{uj), whose contribution to the Liftshitz for- 
mula dominates over any possible magnetic response that 
the metamaterial may have, making attractive a Casimir 
force that would otherwise be predicted to be repulsive. 
We have recently addressed many of these issues in 
the context of the Casimir-Lishitz theory and metamate- 
rials [2l| . The purpose of the present work is to further 
investigate the physics of Casimir interactions between 
metamaterials, focusing on effects not previously consid- 
ered in depth in the Casimir literature, such as optical 
anisotropy in magnetodielectrics and the feasibility of the 
crossover from attractive to repulsive Casimir forces with 
realistic metallic-based and dielectric-based metamateri- 
als. 



their respective polarizations (transverse electric (TE) or 
transverse magnetic (TM)). The operator exp(— /Cd) rep- 
resents one-way propagation between the two plates, and 
has matrix elements 



k|| |e-'='^|k') = cM-dJkl + e/c^) 5^^\h\ K)- (2) 



'II 

When both plates present homogeneity in the plane of 
the interface, only specular reflection takes place, and 
TZj is also diagonal in the transverse momentum basis. 
This means that 



(k|||7^,|^)=R,<5(2)(k||-k'), 



(3) 



where Rj is the 2x2 reflection matrix on the j'-th plate. 
Note that the reflection matrices here are evaluated at 
imaginary frequencies a; = i^, and this requires the well- 
known analytic properties of the permittivities and per- 
meabilities in the complex frequency plane. For general 
anisotropic media these reflection matrices are defined as 



R, 



TE,TE/.^ , N TE.TM/., , x 

TM,TE/-(- , \ TM,TM/.(; , ^ 
^ («C,k||) r- (J^,k||; 



(4) 



II. CASIMIR-LIFSHITZ FORCE BETWEEN 

ANISOTROPIC MAGNETODIELECTRIC 

MATERIALS 

A. The scattering approach 

Techniques for the evaluation of the Casimir force have 
evolved very quickly in the last few years, pa ving the 
way for precise analytical |22l | and numerical [23| cal- 
culations in non-trivial geometries. A particularly ap- 
pealing method is the so-called scattering approach, pi- 
oneered in Casimir physics by Balian and Duplantier 
to compute the free energy of the electromagnetic field 
in regions bounded by material boundaries of arbitrary 
smooth shape [2J|. The free energy is expressed as a 
convergent multiple scattering expansion of ray trajecto- 
ries propagating between the material boundaries. This 
method, first used for perfect conductors, was extended 
to real materials in recent works [25l.l26l|. allowing in prin- 
ciple the computation of the Casimir interaction between 
arbitrarily shaped material scatterers. 

Since a thorough discussion of the scattering approach 
would take us too far afield, we simply present the for- 
mula for the zcro-tcmpcrature Casimir energy per unit 
area A between two parallel plates separated by a vac- 
uum gap of width d: 



E{d) 
A 



fi f ^ log det [1 - 7^le-'='^7^2e-'='^] 
Jo 27r 



(1) 



where TZj = 7?.j(k|[,kJ|,p,p', 



= i^) is the reflection 
operator associated with reflection on the j-th plate 
(j = 1,2). Here k|| and kj, are the transverse wave vec- 
tors {i.e. projected onto the planar interfaces) for inci- 
dent and reflected waves, respectively, and p and p' are 



where r^'^ is the ratio of the amplitudes of a reflected 
field with p'-polarization and an incoming field with p- 
polarization. 

Using Eqs. ([3]) and (|2]) in Eq. ([T]), we get after some 
manipulations 



m ^ „ /-|i 



ri/^>°^*'['-^--^--^'-'i 



and the expression for the force per unit area follows: 



(5) 
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2h 



d^k» 



2Tr J (27r)2 



■ lU Tr 



Ri R2 e 



-2Kad 



l-Ri-R2e-2-^3d' 



(6) 
where K^ = ./k?, +£,'^/c'^. A positive (negative) value of 

the force corresponds to attraction (repulsion). Despite 
the fact that we have assumed homogeneity on each of 
the planar interfaces (which is a reasonable assumption 
when describing metamaterials with an effective medium 
approach), Eq. ([6]) is still fairly general: it may be ap- 
plied to dispersive, dissipative and anisotropic media; all 
that is needed are the appropriate reflection matrices. 

Let us consider the setup depicted in Fig. 1, where we 
have a metallic semi-space occupying the region z < —d 
facing a magnctodielectric semi-space 2 > 0. The reflec- 
tion matrix Ri characterizing the metal-vacuum inter- 
face is given by the standard Fresnel coefficients |27| 



TE.TE 



1 (*'^,k||) 



TM.TM/-^ , \ 



A'3 - ^kf^+e,itOe/c^ 



•Ti («C,k||) =ri ■ (iC,k||) =0, (7) 




z = — d 



z = 



FIG. 1: Typical setup used throughout this paper to compute 
the Casimir-Lifshitz force between a metal and a metamate- 
rial. 



where ei(tj) is the permittivity of the metal. The ele- 
ments of R2 arc only given by Fresnel-like formulas when 
the MM is isotropic, in which case 



TE.TE 



(*C,k||) 



TM.TM 



(*?,k||) 



^^2{iOK3 


V II 


+ M2(JC)e 


.M)e/c^ 


Mi^Ks 


^aAF 


+ M2(Je)e 


.M)e/c^ 


e2mK3 


- y'fc| 


+f^2m^2me/c^ 



TE,TM 



(*?,k||) 



(»e,k||)-0, (8) 



£2 

TM.TE 
"^2 



where 62, ^2 arc respectively the permittivity and the per- 
meability of the metamaterial. 

However, as magnetodielectric MMs can generally be 
optically anisotropic, the reflection matrix R2 for the 
MM- vacuum interface is in general not given by the usual 
Fresnel formulas ([8|). In their most general form meta- 
materials can be bi-anisotropic. rneaning that the consti- 
tutive relations have the form 1281 



D = e E + K H, 
B = C E + ^ H. 



(9) 
(10) 



Here k, and C, arc the magneto-optical permittivities, and 
they describe magnetic-electric cross-coupling. There are 
indeed some metamaterials in which the magneto-optical 
tensors k and (^ are not negligible |23; but since these 
properties can be almost entirely suppressed by using a 
sufficiently symmetric unit cell [30|, we assume hence- 
forth that D = e • E and B = /x • H. We also assume 
again that the material tensors e and // are functions of 
frequency only, neglecting any possible spatial dispersion. 
Even without bi-anisotropy the physics of 
(uni)anisotropic materials is still very rich [3ll . |32| |. 



It is very common to describe them according to their 
degree of symmetry; in crystallographic theory this leads 
to Bravais lattices and their associated point groups 
|33 |. This classification is also very useful for the study 
of metamaterials, since they may usually be described in 
terms of unit cells (split-ring- resonators [3J| , nanopillars 
[13J . nanorods [Ijl, nanospheres [3^, |3al, etc.) arranged 
in a periodic lattice. The most extreme anisotropic 
situation is when the only symmetry of the unit cell is 
inversion with respect to the origin. In this case, known 
as the triclinic system, both the permittivity and the 
permeability tensors have nine non-zero components |32l | 
in a given orthogonal coordinate system, making the 
formulation very cumbersome. Although it is certainly 
possible to diagonalize at least one of the tensors by 
choosing a suitable basis, the angles formed by the 
eigenvectors depend upon frequency in the triclinic 
system [23, [331 ■ Since the force © is an integral over all 
frequencies, this frequency-dependent diagonalization is 
of little help for purposes of calculating Casimir forces. 
Fortunately, it is still possible to investigate anisotropic 
effects in the Casimir force without going into such an 
involved case, so we restrict ourselves in the next two 
subsections to basically two types of anisotropy. 



B. Reflection matrices for uniaxial (out-of-plane) 
planar metamaterials 

In this subsection we calculate the reflection matrix 
R2 for the case of a planar interface between vacuum 
and a uniaxial magnetodielectric medium that is isotropic 
on the interface plane, i.e., whose electric and magnetic 
anisotropic directions coincide and are perpendicular to 
the interface. In optical terminology, this is an example 
of a uniaxial medium [27| with the optic axis coinciding 
with the anisotropic direction. It is known that for uniax- 
ial lattices, that is, the ones belonging either to the trig- 
onal, tetragonal and hexagonal crystallographic systems 
[33J . the electromagnetic tensors are diagonal in the co- 
ordinate system defined by any two orthogonal directions 
in the symmetry plane and the optic axis [37[. Therefore, 
choosing the interface as the xy plane and the anisotropic 
medium to be the half-space defined by 2 > (see Fig. 
2), the permittivity and permeability tensors are given 
by 



e. 








^J^.xx 











; ^i^j = 





l^xx 





zz 










tJ^zz 



,(11) 



where we have used Cyy = e^x and fiyy = fixx-, whose fre- 
quency dependence is implicit. Although the calculation 
of the reflection matrix for a metamaterial with a single 
out-of-plane anisotropic direction is relatively simple and 
akin to the isotropic case, for the sake of completeness 
we briefly review this calculation, which is relevant to 
several metamaterials having such anisotropy [38|. 




z=0 



FIG. 2: An incident plane wave impinging on a uniaxial meta- 
material with its optic axis perpendicular to the 2 = plane. 



Let us assume that a plane wave with wave vector 
k and polarization p impinges upon the interface from 
z < Q (region 3, vacuum) towards the metamaterial (re- 
gion 2). Given the rotational symmetry about z^ without 
loss of generality we can choose our coordinate system so 
that the plane of incidence (defined by k and z) coincides 
with the xz plane (see Fig. 2). In order to solve the 
reflection-refraction problem we have to know how waves 
propagate in the anisotropic medium. In this particular 
case of uniaxial anisotropy orthogonal to the interface it 
may be shown by direct substitution that TE waves 



E 



TE 



Eaye 



i{kj^x+k^z) —iujt 



'^2,x , '^2,z 


u? 


— — H = 


= C2 '^^ 


Mzz ^J■xx 



are solutions to Maxwell's equations provided that 

[in the MM] 

A:? „ + A:? , = ^7 [in vacuum] . 

In a similar fashion, TM waves 

H™ = ifoye^<'=-"+'=^")e-*'^*, 
are solutions to Maxwell's equations provided that 



'^'l,x 



'^2,z 



-fj.xx [in the MM] 



A-2 -I- A-2 



in vacuum 



(12) 

(13) 
(14) 

(15) 

(16) 

(17) 



Therefore there is no polarization-mixing, and conse- 
quently the off-diagonal elements of the reflection matrix 
vanish: rj ' (i^,k||) = r2 ' (j^,k||) = 0. This al- 
lows one to consider separately the reflection of TE and 
TM waves. 



Let Ein = ^oye*('=i>-^+'=i--^) be a TE field incident 
from the vacuum side. Given the translational invari- 
ancc of the material properties along the planar in- 
terface, only specular reflection occurs, which implies 
that both X and y components of the wave vector k 
are continuous. Therefore, the reflected TE field is 
Eref = r2£;oye*(''i-'="''i-'^\ and the transmitted TE field 
is Et = t2Soye'<'=^'-"+''^-"', with q ^ = {uo^l^)e^,iixx - 



kf ^fJ.xx/ l^zz- Imposing the bomidary conditions on the 
TE modes, wc have 



E, 



in,J/ 



Ei-cty 



Et.v ^ l+r2=t2 



Hi, 



Hrci,x = Ht^x => (-1 + ?'2)/si,2 = -t2- 



k2,z 

'. : 

^xx 



from which it follows that r2 = {fJ-xxki^z — k2,z)/itJ-xxki^z + 
k2.z)- Evaluating this expression along imaginary fre- 
quencies Lo = i^, one obtains the TE-TE reflection am- 
plitude on the vacuum-MM interface [33 : 



^2„„i «>k||) 



fJ'XxKs 



f^xx u2 _]_ , , p s" 
u- II '^ i^xx'^xx'^ 



li-xxK-i -\- 



i^xx u2 [ II c s_ 

A,|| ^^ I^XX'^XX^ 

fJ'ZZ II L, 



where, we recall, /vs 



H 



$2/c2 and fcf? = kl 



(18) 



ky. 



Following similar steps, the TM reflection amplitude on 
the vacuum-MM interface can also be derived: 



TM,TM/.^ , N 



nA', 



^XX htZ 



^Xxf^XX ^2 



rK. 






(19) 



^xxf^xx "TH 



C. Reflection matrices for biaxial, anisotropic 
magnetodielectrics 

In ascending order of symmetry, the crystals belonging 
to the triclinic, monoclinic and orthorhombic crystallo- 
graphic systems [33| are known as biaxial crystals, since 
they are characterized by two optic axes. In this subsec- 
tion we shall restrict ourselves to the orthorhombic case 
j40| , which allows simultaneous diagonalization of e and 
/x in an orthonormal basis. The calculation of the reflec- 
tion matrices for the other two types of biaxial metama- 
terials is conceptually equivalent but more cumbersome 
since the material tensors cannot be brought to diagonal 
form in a frequency- independent basis. 

Let us then consider the system described in Fig. 3, 
which is similar to Fig. 2 but with an orthorhombic meta- 
material on the right side. Assuming it is possible to pre- 
pare the MM in such a way that one of the eigenvectors 
is perpendicular to the interface, then the diagonal basis 
is just {x, y, z} and the electromagnetic tensors are given 
by 



(20) 



<^xx 










[i-xx 











^yy 





; M^i = 





Ihv 











^zz 










i^zz 




z = 



FIG. 3: An incident plane wave impinging on a biaxial meta- 
material with orthorhombic symmetry (see text). 



Several mctamatcrials can be described by material ten- 
sors like (PP]) : a good example is the fishnet design used 
in M- 



Mctamatcrials with two optic axes, even those with 
the simplest orthorhombic symmetry, are much harder 
to treat than those with an out-of-plane, uniaxial optic 
axis described in the previous subsection. The reason 
is that Maxwell's equations do not support transverse 
waves for biaxial materials: neither TE nor TM waves 
are solutions inside the material, and the off-diagonal el- 
ements of the reflection matrix do not vanish. This also 
happens for uniaxial materials with in-plane optic axes, 
whose reflection matrix can be obtained as a particu- 
lar case of orthorhombic materials with Syy = e^^ a-nd 
fiyy = fizz- The Casimir interaction between two dielec- 
tric semi-spaces with one in-plane optic axis was treated 
in [41| and used in the experimental proposal to measure 
the Casimir torque between birefrigent plates [42 |. 

The calculation of the plane-wave solutions to 
Maxwell's equations is simplifled using a coordinate sys- 
tem attached to an incident wave from the vacuum side. 
Let a plane wave with incident wave vector k impinge 
on the interface forming an angle 0^ with the normal 
direction (see Fig. 3). Let {x',y',z') be the coordinate 
system attached to the corresponding plane of incidence, 
that forms an angle ip with the x axis. The optical tensors 
in this new coordinate system are 



J 



£i'i 



Cxx cos^ ip + eyy sin^ ip (e^x - ^yy ) sin ipcosip 
{exx ~ f^yy) simp cos ip exx sin^ ^ + eyy cos"^ ip 











and 



fJ-i'j' 



flxx COS^ ip + fiyy Sin^ ip {fixx - fJ-yy) simp COS ip 
ifj'xx - fJ-yy) sin ip cos ip jlxx sili^ ip + jJLyy COS^ ip 

fizz 



r 



The expressions for the incident fields are 



Ei, 



H;,, 






e™y'-e™-(<Zi„i'-fc.-z') 



-,i{k^i x' +qi^z' —Lot) 

(21) 



(22) 

where e™, e™ are given amplitudes and we defined 
kx' = (w/c)sin0in and gin = (w/c) cos6'in. The reflected 
wave has a similar expression: 



E, 



Hic 



ejjy' 


-e™-(gi„i' + fc,,z')" 

u; J 


i{k^^/ x' —QinZ' —cut) 

(23) 


'el^if 


+ e;^J-(ginx' + /.vz') 


gi(fc^,a;'-9i„z'-wt) 
(24) 



where we have used q^ci = — ffin- Our problem now con- 
sists in flnding the amplitudes e™ e™, so we can con- 
struct the reflection matrix (|4]). In order to obtain the 
reflection amplitudes, however, it is necessary flnd the 
transmitted fields as well, which means that we have to 
solve Maxwell's equations in the metamaterial. 
Let us assume plane waves 

E = e(z')e^('=^'^'""*) ; e = (e,., e^., e,.), 
H = h(z')e^^'=='"'""*^ ; h={hx',hy,,hz,), (25) 

as solutions to Maxwell's equations in medium 2, where 
we have already deduced the x' dependence from the 
phase-matching condition on the interface {kx' is con- 
served across the interface). By substituting ([25|) into 
the Faraday and Ampere-Maxwell laws 



V X E 






V X (/2" 



B) 



1 dje ■ E) 
c dt 



(26) 



respectively, we see that the z' components can be chm- 
inated as 

e^/ = -ck^ihyijuje^i^i ; h^j = ck,j,iey> jujp.:,! ^i . (27) 

In order to determine the remaining x' and y' components 
of e and h it is convenient to introduce a vector u with 
components u\ = e^' , U2 = Sy', u^ ~ h^' and U4 = hyi. 
With the ansatz Uj = Uj(0)e**^ we obtain the foUowing 
hnear system of equations: 



L • u = q u, 

UJ 

where the non-zero elements of the matrix L are: 
Lid. = -L2i = -(Mxk - ^iyy)sm.ipcoii^p, 



(28) 



^14 



2 ^Jiyy cos^ ip - ^ixx sin^ (^, 



L23 = tJ-XX COS^ ip + ^iyy SIT? if, 

L31 = -L42 = {e^x - f-yy) sin ip COS LP, 
c 



-^32 — 



^^ 1- £.„.„ cos^ ip + ej;j, sin^ ip, 

UJ^lJ.zz 

L41 ~ ^€xx COS (p> — e^j, sin Lp. 



(29) 



The condition for non-trivial solutions (det(L + Luq/c) = 
0) gives us the equation that determines the possible val- 
ues of (7, namely 



^,^.a)(^,^-b)^c. 



(30) 



m = 1, 2, we may write the general solution for u as 

u(z') = J2 u(™)(0)e'''"'"' 



m=l,2 

E 

m=3,4 



U 



(m) 



(0)e- 



(32) 



where we have used q^^^ = — g*^^^ and q^'^'' = ~q^^\ It 
is easy to see that the refraction of a wave coming from 
z' < can only give rise to positive propagating waves 
(in the sense defined above), from which we conclude that 
u(3) (0) = u(4) (0) = 0. Therefore, the transmitted field 
into the anisotropic magnetodielectric medium is 



Et 
Ht 



*(fex 



m=l,2 



(0) e*« 



(33) 



In order to find the amplitudes u'"'^ (0) we have to im- 
pose the proper boundary conditions on the fields. In this 
case they just require the continuity of E^i , Eyi , H^' , Hyi 
across the interface. Using (PT1) - ([M)) and (|55)). one derives 
the following boundary conditions: 



9in(e 



TM 



-gin(e 



„TE 



TE 



„TM 



where 



TMn 
^rcf ) 



„TE 



°ref I 



TM 



m=l,2 



ei"^(0), 



m=l,2 

^ E e^(o)' 

m=l,2 



(34) 



A = Li^Lsi + L14L41, 

B — L23L32 + L24L42, 

Ci = L13L32 + L14L42, 

C2 = 2723-^31 + -^24^41, 

C = C1C2, 
whose four solutions q^^'^^ [m ~ 1, 2, 3, 4) are 



w 1 



qira) ^^"^ Ja + B±^{A- Bf + AC. (31) 
c V2 

These solutions may be conveniently split into two pairs, 
according to the sign of Re q^™) - solutions with Re 
qr("i) > Q (j^e (7^™) < 0) define positive (negative) prop- 
agating waves. If we denote the positive solutions by 



This system of equations is unsolvable as it stands, given 
the large number of unknowns. It is possible, however, 
to use ((28|) to express all the transmitted amplitudes in 



terms of just one, say e^7 : 



v(™) 



^(™) 



v(™) 



(m 



^? 



(0) (g(™))2 - (u;7c2)A 



(0) 
(0) 



(^Vc2)Ci 



UJ L31 



(0) 
(0) 



cq 



(m) 



--^a(™), 



cq 



(m)' 



w L41 UJ L 



e'j\Q) 



c q 



(rn) 



cq- 



(m) 



Using these definitions, we can rewrite 



as 



-1 a(i) «(2) 

cgi„/cj -/3(i) -/3(2) 

C(7in/cJ 1 1 

-1 7(1) 7(2) 



„TE 
^ref 
pTM 

4V(o) 

4^^(0) 



/ TM 



(35) 



In order to find the refleetion eoefRcients, we must solve 
(p5|) for tlie reflected amplitudes. For the sake of clarity, 
let us do this separately for e™ = 0, e™ 7^ and for 
e™ ^ 0, e™ = 0. In the first case, Cramer's rule 
immediately yields 
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and in the second case we have 
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where M is the 4x4 matrix in (|55|) and 
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III. METALLIC-BASED METAMATERIALS 
AND THE CASIMIR EFFECT 

Metamaterials may be roughly divided into two classes. 
The first class consists of MMs that are partially or to- 
tally based on metallic structures. In this section we con- 
centrate on these metallic-based MMs, which were pre- 
viously considered by us in [2l|. We study in detail the 
effects of optical anisotropy on the Casiniir-Lifshitz inter- 
action with magnetodielectric media. The second class 
consists of MMs based purely on dielectric materials, that 
we shall treat in the next section. 



A. Isotropic metamaterials 

Before going straight to the calculations, it is necessary 
to point out that metallic MMs may be also divided into 



two types, which we shall characterize as (i) connected 
and (ii) non-connected. As the name suggests, in the 
connected MMs the metallic part is partially or totally 
interconnected throughout the metamaterial [l5[, while 
in the non-connected it is not [ij, [lj| . This distinction 
is important because in connected MMs there is a net 
conductivity contribution to the dielectric function due 
to the metallic part, while in the non-connected MMs the 
background is effectively non-conducting. 

Let us begin with the simple example of a metallic 
half-space 1 in front of an isotropic, connected metallic- 
based metamaterial 2. For the metal we assume the usual 
Drude model 
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(41) 



where Vli is its plasma frequency and 71 the dissipation 
coefficient. For the second half-space we have to be more 
specific about the MM we want to consider. In the sim- 
plest description isotropic, connected metallic metamate- 
rials may be described by a dielectric response accounting 
for both a resonance and a Drude contribution: 
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(42) 



where fie, coe and je are respectively the effective electric 
oscillating strength, the resonance frequency, and the ef- 
fective dissipation parameter of the resonant part, and 
VId and jd are the Drude parameters of the metallic 
background of the MM. The filling factor / roughly quan- 
tifies the fraction of metallic structure present in the MM. 
The magnetic permeability is given by a resonant part 
alone: 



fJ-2iuj) = 1 - 
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where J7,„, Wm and 7^ are defined analogously to their 
electric counterparts. In Fig. 4 we plot the Casiniir- 
Lifshitz force between a metallic half-space and an 
isotropic metallic-based planar metamaterial described 
by ((42|) and ((43|) for different filling factors at zero tem- 
perature. We see that without the Drude contribution 
(/ = 0) there is repulsion for a certain range of dis- 
tances, as long as the half-space 2 is mainly magnetic 
(£2(*C) < M2(*0)- However, as we "turn on" a metal- 
lic background (/ > 0), the permittivity grows stronger 
and reverts the previous relation for a larger and larger 
range of frequencies, up to the point where the magnetic 
activity is no longer able to produce repulsion. 

An idealization carried throughout the paper is that 
both the metal and the metamaterial are infinitely long 
in the z-direction. When we have slabs of finite thickness 
instead of half-spaces, the reflection coefficients change to 
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where Kj — ./ku + ^j{i^)ej{i^)^'^/c'^, dj is the thickness 

of the j-th slab, and we are assuming that both slabs are 
surrounded by vacuum. From the previous expression 
we see that corrections to the half-space reflection coeffi- 
cients (at imaginary frequencies) are exponentially small 
when both products Kidi and K2d2 are sufficiently large. 
This basically tells us that estimates on lower bounds for 
di and d2 are actually model dependent (given that Kj 
depends on the properties of medium j), so in order to 
discuss those estimates we have to be more specific. For 
a metal described by (|4T|) and a wave arriving at normal 
incidence, we have 
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For 



vp — 27rc/r2i, Xd = 2TTc/"fi and A = 27rc/^. 
high frequencies, we have A <C Ad and then p5|) becomes 
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For typical metals Ap/27r is around 10 — 20nm, so, at least 
for high frequencies, the contribution to the Casimir force 
of a slab some tens of nanometers wide approaches the 
contribution of a half-space. However, in the opposite 
limit (A » Ad), we have 
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FIG. 4: The ratio F/Fq for a gold half-space facing an 
isotropic, interconnected and silver-based metamaterial. F/A 
is the Casimir force per unit area in this setup, Fc/A = 
hc-K^ /2AQa*' is the Casimir force per unit area between two 
perfect plane conductors, and _F < 0(-F' > 0) corresponds to a 
repulsive (attractive) force. The frequency scale Q. = 2-nc/ A. 
is chosen as the silver plasma frequency Q,d = 1-37 x 10^® 
rad/sec. Parameters are: for the metal, f7i/Sl = 0.96, 
7i/r2 = 0.004, and for the metamaterial, fio/fi = 1, 7D/fi = 

0.006, Q.,JQ. = 0.04, ^m/0. = 0.1, LOe/0. = Lj,„/n = 0.1, 

7e/fi = 7m/f^ = 0.005. The inset shows the magnetic perme- 
ability ^2 (iS,) and the electric permittivity £2 (i£,) of the MM 
for the different filling factors. 
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where ctq = ^1/7 is the static conductivity and 5{^) is 
the skin depth of the metal at imaginary frequencies. 
Thus, for long wavelengths we see that d scales with -s/A, 
leading to the conclusion that the half-space approxima- 
tion is not good for sufficiently low frequencies. Fortu- 
nately, for typical materials this is no source of concern, 
since the integration range where A 3> A^ holds is very 
small compared to the effective integration range, allow- 
ing us to push the slab approximation up to very small 
frequencies with almost no effect in the final result. One 
might wonder what happens for oblique incidence, but it 
is easy to see that the more oblique the incident angle is 
the better the estimate for di holds, since Ki gets larger 
and larger (this only means that reflection gets easier as 
the incidence angle gets larger, as physically expected). 
The effect of finite thickness in the Casimir effect was the 
object of several papers [43, notably in [4^ where a sys- 
tematic procedure was developed to deal with any given 
number of arbitrary slabs. The effect of finite thickness 
was also studied in the specific context of Casimir force 
and metamaterials [1^, [13, S^] , where it was found that 
having a layer of a MM instead of a half-space reduces 
the intensity of the repulsion force and also the range of 
distances where it occurs. 



Uniaxial Metamaterials 



Electric anisotropic effects in the Casimir interaction 
have been thoroughly studied in the literature [41, 42, 43, 
|49| . l5g | , but until recently there was no compelling reason 
to study the consequences of magnetic anisotropy. This 
changed with the advent of metamaterials, and an in- 
vestigation of magnetic anisotropy is now in order. The 
best place to start is to consider uniaxial out-of plane 
metamaterials, since they constitute the simplest depar- 
ture from the isotropic case. This type of anisotropy is 
quite common since it arises naturally when a material 
is built as a stack of different layers, as is the case for 
several kinds of MMs [la, [3a, [mI, [SJI . We are particu- 
larly interested in the case where the resulting medium is 
characterized by different degrees of conductivity in the 
plane of symmetry and in the perpendicular direction to 
it. 



Let us begin by characterizing the electric and mag- 
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FIG. 5: The effects of uniaxial anisotropy in the Casimir force 
between a gold semi-space and a metallic-based connected 
MM with weak Drude background. The distance is fixed to 
d = A and repulsion corresponds to negative values of F/Fc ■ 
All parameters are the same as in Fig. 4 except for the filling 
factors fx and fz, which are the variables in this plot. 



netic properties of our uniaxial metamaterial: 
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(48) 



where the different filling factors fx and fz account for 
the possible anisotropy in the metallic character of the 
MM. As we have seen in subsection II-B, the reflection 
matrix for such a metamaterial is diagonal and given by 
([T8| and ([T9|) . For metallic-based mctamaterials with 
large in-plane electric response exxi^S.) 3> 1 at low fre- 
quencies, it is clear from Eqs. (UHl ^^ that anisotropy 
plays a negligible role in the determination of the reflec- 
tion coefficients when there is a dominant Drude back- 
ground. In order to better appreciate the effects of 
anisotropy we assume henceforth a small or vanishing 



Drude contribution. In Fig. 5 we show the Casimir 
force for a metamaterial that has only electric anisotropy 
ifJ-xx = l^zz), which is completely coded in different filling 
factors {fx 7^ fz, afi other parameters being the same). 
We see that a repulsive force (F/Fq) arises only for con- 
siderably small values of both fx and fz, from which 
we conclude that killing the Drude background in the 
z-dircction alone is not enough to produce Casimir re- 
pulsion. 



C. Biaxial metamaterials 

Continuing our track towards more complicated me- 
dia, we now tackle the biaxial orthorhombic case. Let us 
consider a metamaterial characterized by the following 
dielectric and magnetic functions in the basis defined by 
its eigenvectors (see subsection II-C): 
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We are particularly interested in the case where Sxx is 
close to eyy but in general significantly different from 
ezz- This means basically that the MM is only slightly 
anisotropic in the plane of incidence. Our motivation 
in studying this particular limiting case is that it is a 
good approximation for certain types of metamaterials, 
such as those based on fishnet designs [l^ . Note that we 
are already assuming magnetic in-plane isotropy, which 
is consistent with a small electric in-plane anisotropy. We 
may then rewrite the material tensors as 
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(50) 



where S{uj) = {eyy{uj) — exx{^))/exx{^) <C 1, and perform 
the calculations only up to first order in 5. The evaluation 
of the determinants in ([5S)) - ([551) requires the knowledge of 
matrix elements of L defined by equation (j29|l and also of 
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the solutions of equation (pO|l . This last step is simplified 
at first order in 6 because C/e^^. ^ {e^x ~ ^yyY l^tx ~ ^^ ■ 
Therefore 



c c 



(51) 



and then the a^™) coefficients reduce to a^^^ = and 
q;(2) = (i? — A)/Ci. Performing now some straightfor- 
ward calculations we get, up to 0{5^), 

TE.TE/ N TE.TE/ \ , e TE,TE/ \ 

^2 ('^)=^2„„i (^) + '5^2.1 ('^). 



TM,TE 



TE.TM 



TM,TM 






Sexx^J.xxqtlnqin sin 2i^ 



(gtm + gto)(gtm + £xxqin){qtc + fJ-xxqm) ' 

SexxfJ-xxq-m sm2(p 

(qtm + ^xxqin)iqtc + y-xxq-m) ' 



TM,TM 



TM,TM/ 



/ \ I C liVi,±iVi/ \ 



(52) 



where 
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are given respectively by ([TS]), 



T9|l . and we have also defined 
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FIG. 6: The Casimir force between a gold half-space and an 
orthorhombic, slightly in-plane anisotropic MM for different 
values of the filling factors f^ and fy . The bands are charac- 
terized by a certain value of fx, as shown in the legend, and a 
continuum of values of fy, from fy — 0.8 fx to fy = 1.2 fx- All 
the other parameters involved are exactly the same as those 
used in Fig. 4. 
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Now let us return to the general structure of the Lif- 
shitz formula. Since R2 is not diagonal we should expect 
contributions coming from the non-diagonal terms in ([5|), 
but it can be shown that they are all 0(i5^), and therefore 
can be dropped. After a few rearrangements we arrive at 
our final expression for the Casimir pressure: 
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An easy consistency check of this result is to take the 
zero-anisotropy limit, which reduces immediately to the 
uniaxial result, as it should. A less trivial result is to 
obtain the non-retarded limit of expression ([53]) . which 
can be shown to be consistent at first order with other 
results in the literature fill. IJtI. list . 

Let us now assume that all the in-plane anisotropy 
is coded in the filling factors, just like the out-of-plane 



anisotropy was in the uniaxial case. In Fig. 6 we show 
the effects of a slight in-plane anisotropy on the Casimir 
force. Each band in the plot corresponds to a different 
value of fx, and its width is given by a ±20% variation 
of fy around fx- We see that the anisotropy effect is 
more pronounced at small distances, because in the non- 
retarded limit the contribution of the electric response to 
the Casimir force is maximized. 



D. Dissipation effects 

Let us now turn from considerations of anisotropy to 
other practical issues for MMs and Casimir interactions. 
It is known that dissipation plays an important role in 
metallic-based metamaterials, especially those operating 
at high frequencies. In Fig. [7| we show the effect of a 
sinrultaneous modification in the electric and magnetic 
dissipation coefficients; it may be clearly seen that an 
equal change in the rates 7e/i^e ^^nd "fm/^m favors at- 
traction. In the insets (a) and (b) we show respectively 
the effects of changing only the electric and magnetic 
dissipation, that may be straightforwardly interpreted in 
light of the discussion presented in [iO|. From ((42|) we 
see that an increase in je makes e smaller, pushing the 
metamaterial slightly closer to the Boyer limit (that is, 
£1 -^ 00, ^1 = 1, £2 = 1, ^J■2 -^ oo)- We should thus ex- 
pect an increase in the Casimir repulsion as we make 7e 
larger, and that is exactly what is observed in the inset 
(a) . A similar reasoning nray be applied to inset (b) , but 
since this time we are going away from the Boyer limit, 
repulsion diminishes as we increase 7^- 
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FIG. 7: The ratio F/Fc between a gold half-space and an 
isotropic silver-based metamaterial for different values of the 
dissipation parameters. The main plot shows the effect of the 
simultaneous modification of electric and magnetic dissipa- 
tion. Inset (a) shows the effect of electric dissipation alone for 
different values of the ratio ye/^e— 0.1 (solid), 0.5 (dashed), 
2.5 (dotted). Inset (b) shows the effect of magnetic dissipa- 
tion alone for different values of the ratio ym/^m= 0.1 (solid), 
0.5 (dashed), 2.5 (dotted). The filling factor is / = 10"* in aU 
three plots, and all other parameters except the dissipation 
coefficients are the same as in Fig. 4. 



E. Temperature effects 



isotropic materials at a finite temperature T 
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where the prime in the summation means that the 
n = term is multiphed by 1/2, /3 = l/fceT, K3 = 

Jk?. + Cn/c^i £,n ~ 2T:n/hp are the Matsubara frequen- 
cies, and the reflection coefficients are given by ([7]) and 
([8]) with ^,1 instead of ^. From (|55p we see that for large 
distances (provided that fc^Td/Zic ^ 1) the n = domi- 
nates all the others, and wc may approximate the Casimir 
force by 
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where k = |k||| and the reflection coefficients are evalu- 
ated at the zeroth Matsubara frequency Co = 0- The key 
difference from the setups using dieletrics or Drude met- 
als to the one with plasma metals is that in the former 
cases we have lim^^o e(C)C^/c^ = 0, leading to 
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k-k 

k + k 
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while in the latter we have 



The effects of temperature in the Casimir force be- 
tween metamaterials and dielectrics were thoroughly dis- 
cussed in [ifll, where the authors show that in this case 
temperature works against repulsion, or in other words, 
that for sufficiently high temperatures repulsion is com- 
pletely overturned into attraction. In this section, we 
wish to extend that discussion for a metallic plate fac- 
ing a MM and compare the situations where the metal is 
modeled by either Drude or plasma permittivities. 

In Fig. 8a we show the Casimir force for different tem- 
peratures between a Drude metal and an isotropic MM 
with no Drude background. In this case we see that tem- 
perature also works against repulsion, but in such a way 
that keeps the repulsion window quite open for tempera- 
tures as high as r = 600K, allowing for repulsion at room 
temperature, at least in principle. Something even more 
interesting happens when we change the Drude metal 
by a plasma metal (i.e., vanishing relaxation parameter 
71 = in eq. (J4T|) ). as shown in Fig. 8b. In this case, we 
see that not only a temperature increase does not switch 
back the force into attraction for large distances, but it 
actually increases repulsion in that regime. It is possible 
to explain this phenomenon in simple terms using the 
Lifshitz formula. Let us consider the force between two 



7-1 (o,fc) 



yP 



k 4- ^fc2 + ni 



1 <0. 



(58) 



This means that for dielectrics or Drude metals facing 
a MM, the only contribution to ([56|) comes from the 
TM zero mode, which is always positive (given that 
r-^ {\),k)r2 \S>,k) > 0). Smce this term domi- 

nates for large distances, we conclude then that the force 
is attractive in this regime. However, for plasma metals 
facing a MM we see that both TE and TM zero modes 
contribute, and while r^ ' (0, k)r2 ' (0, k) is posi- 



tive the product r^ ' (0, k)r2 ' (0, k) is not, due to 
the different signs of r^ ' (0, fc) < and of rj ' (0, k) 



^J■2(0) + 1 



(59) 



We see then that the sign of the force depends on a deli- 
cate balance between the TE and TM contributions, and 
it so happens that for our chosen parameters the TE term 
overwhelms the TM term and repulsion is sustained for 
all distances above the crossover from attraction. The 
fact that repulsion is enhanced is also easily explained, 
since a simple analysis shows that for large distances (|56p 
may be put in the form C/f3d^, where C is a constant 
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FIG. 8: Temperature dependence of the Casimir force be- 
tween a metallic plate and a metamaterial. We plot the 
Casimir force between a metamaterial and a Drude metal (a) 
or a plasma metal (b) for different temperatures. We stress 
that negative values of the force characterize repulsion, and 
that all parameters are the same as the ones used in the / = 
curve of Fig. 4. 



the currents in each pillar are running in opposite direc- 
tions) at a certain frequency, where the electric dipole 
contributions of both pillars are nearly canceled out and 
the effective current loops produced by the pairs give 
rise to magnetic dipole contributions, resulting in a non- 
trivial magnetic activity. 

Unfortunately, a detailed treatment of the metama- 
terial previously described is beyond the scope of the 
present paper. It is still possible, however, to capture 
some effects of metallic non-connectedness and geomet- 
rically built-in resonances through the use of an appro- 
priate toy model. In order to address the first issue, we 
consider a simple MM model consisting of identical, small 
metallic spheres of radius a regularly distributed in a 
host dielectric (non-magnetic) medium. Assuming that 
the metal and the dielectric are characterized respectively 
by the permittivities 
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and that the metallic spheres can be considered in a first 
approximation as electric dipolcs, one can connect the 
medium effective permittivity e2,nc(i^) to the electric po- 
larizability a(a') of a given sphere through the Clausius- 
Mossotti formula [491 
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where / is the metallic filling factor and we have sup- 
pressed the w-depcndence for simplicity. It is also pos- 
sible to show that when {uj/c)a <C 1 (the spheres are 
much smaller than the radiation wavelength), the elec- 
tric polarizability may be given in terms of the dielectric 
function of the metal by the similar relation 



depending on the materials used. It is clear then that if 
C is negative, a temperature increase can only enhance 
repulsion. Our findings for temperature effects using ei- 
ther Drude or plasma models for the metallic plate are 
consistent with the conclusions of Hal . 



a __ £2, mot — e2,d 

a? e2,mot + 2e2,d ' 
and by eliminating a in (j61|) we get 
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This result is known as the Maxwell Garnctt approxima- 



F. Metamaterials based on isolated metallic 
structures 

There are several examples of metamaterials where the 
metallic part is distributed in a non-connected way. In 
jlSt . for instance, the authors put forward a MM consist- 
ing of pairs of metallic nanopillars, regularly distributed 
on top of a dielectric substrate. The pairing of pillars 
is necessary to create an antisymmetric resonance (when 



after the physicist who 
A brief analysis shows 



tion for the permittivity |35l . 146 
derived it in the early 1900s [54 1. 
the main effect of having isolated metallic pieces: the pre- 
vious formula tends to a finite value in the zero frequency 
limit, unlike ([42]) . which describes a connected metallic 
MM. The effective permeability can be dealt in a similar 
way, and it is possible to show that in this approximation 
we have simply ^2,nc(w) = 1. 

As noted earlier, formula ([SH]) accounts only for effects 
of metallic non-connectedness. In order to include the 



13 




FIG. 9: The ratio F/Fc for a gold half-space facing a 
isotropic, non-connected and gold-based metamaterial. The 
parameters for the metal are f22,met/f2 = 0.96, 72,met/f2 = 
0.004, and for the metamaterial we have fie/n = 0.34, 
n^/n = 0.064, LJc/fi = 0.2, Lu^/Q, = 0.15, 7c/n = 0.04, 
7m/n = 0.02, / = 0.1. The inset shows the permittivity and 
permeability inside the MM, as given by (|64p. but as functions 
of imaginary frequencies 5. 



built-in electric and magnetic resonances [55|, we sim- 
ply assume their existence in an ad hoc manner and add 
their contribution to e2,nc('^) and H2,nc{i^), respectively. 
Assuming those resonances can be modeled by Drude- 
Lorentz formulas, we have, finally, 



e2(w) 


= 


e2,b 


(w) + e2,rcs 


{LO) 






= 


(1 + 2/)e2,„,et 


+ 2(1- 


- /)e2,d 




(1 


— /)e2,mot 


+ (2 + 


/)e2,d 




+ 




^l 








^2- 


- LU"^ + ijcUJ 




//2H 


= 


1 - 


^l 






UJ^ - ^2 -1 


- i7mW 





(64) 



The results for the Casimir force are shown in Fig. 9. The 
parameters for the resonant parts e2,res(w) and /X2,rcs('jj) 
are roughly based on the experimental results given in 
[131 for a MM consisting of metallic nanopillars covered 
with a thin layer of glycerine. As indicated earlier, our 
intention here is not to provide a precise description of 
such experiments, but only to estimate how this type of 
metamaterial affects the Casimir force. The embedding 
dielectric, glass BK7, is quite well described by (|60p with 
the parameters N = 3, i^2 1/^ = 1-84, uj2 1/^ = 1.81, 
172,2/^ = 0.47, 0J2.2/^ = 0.28, ^2.3/^^ = ^2,3/^ = 0.014, 
72.1/^ — 72,2/^ — 72,3/^ = 0. It is clearly seen that no 
repulsion is achieved, and the reason is that the magnetic 
resonance created by the MM geometry is too weak to 
overwhelm the electric background. In other words, the 
MM is mainly dielectric, leading to an attractive force. 



IV. DIELECTRIC-BASED METAMATERIALS 
AND THE CASIMIR EFFECT 

Metamaterials based exclusively on dielectrics [3^, [5^, 
|57| are an interesting alternative to metallic MMs. For 
one thing, they provide new possibilities for the construc- 
tion of negative index materials [58| , since they allow for 
both the permittivity and the permeability to assume 
negative values in bandwidths that may be out of reach 
with metallic-based MMs. In addition, dielectric-based 
MMs might be interesting for Casimir force studies for 
the same reason that non-connected metallic MMs might 
be: they do not present a Drude background at low fre- 
quencies, and this is advantageous for the observation of 
magnetic effects in the Casimir force. 

The dielectrics most commonly used in the construc- 
tion of MMs are "polaritonic" crystals [5^ characterized 
by the dielectric function 



epoi(w) = 



pol 



^pol 



+ ^pol+*7polW^ 



(65) 



where Wpoi is a characteristic resonance of the system, 
£00 is the permittivity at very high frequencies, and 
^poi = '^poi-\/e(0)/eoo. In order to fix ideas, let us 
consider a MM made of a regular array of polaritonic 
nanospheres of radius a embedded in an isotropic dielec- 
tric and non-magnetic host characterized by a dielectric 
function e^ ■ For sufficiently long wavelengths and sparse 
arrays, meaning x = ujR/c ^ 1, it is possible to use the 
so-called extended Maxwell-Garnett theory [35| to eval- 
uate the dielectric and magnetic properties of the meta- 
material, giving |35l.[36l| 



t^Ging 



(w) = e/i- 



3ifai 



ifai 



fJ-emgi^) 



x^ - 3ifbi 



(66) 



where / is the array filling factor and ai, 61 are respec- 
tively the electric and magnetic dipole coefficients of the 
scattering matrix of a single sphere, given by [23, l60| 

__ ji(a:poi)[a:ji(a:)]'£poi - ji{x)[xpaiji{xpoi)]'eh 
hy' (x) [xpoiji (a;poi)]'e,i - ji {xpoi)[xhY' (x)]'epoi 

jl {Xpo\)[xjl (x)]' - ji (x) [Xpoljl (2;pol)]' 



bi 



,(+)/ 



hY^ {x)[xpo\jiixpo\)]' - Jiixpo\)[xh\ '{x)]' 



(67) 



where ji{h^) is the spherical Bessel function (Hankel 
function of the first kind) of order one, Xpoi ~ ^Cpoix 
and the prime has the usual meaning of a derivative with 
respect to the function argument. The important thing 
to notice here is the fact that /icmg niay present several 
resonances even when the nanospheres are purely dielec- 
tric, from which we conclude that in this framework we 
do not have to assume an ad hoc resonant behavior; it is 
already built into the theory. 
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FIG. 10: The permittivity etruo(*C) ^'^id permeability /itruo(JC)- 
The parameters are Coo = 2, ^^oij^ = 0.4, tOpoi/^ — 0.15, 

-fpal/fl = 0.001. 



The usual procedure at this point would be to rotate 
expressions (j66p to the imaginary frequency axis and sub- 
stitute them into the appropriate reflections coefflcients, 
but in this case things are not so straightforward. Im- 
plicit in the Lifshitz formula for imaginary frequencies it 
is the assumption of analyticity of e{u;) and /i(w) in the 
upper half-plane, a condition that eemg('^) and /icmg('^) 
do not satisfy. In order to overcome this obstacle we have 
to remind ourselves that expressions (|66[) were derived as 
approximations to the true permittivity e(Lo) and perme- 
ability /.t(w) only for a given range of real frequencies, 
namely, for tu such as u>R/c <C 1 . This means that while 
e{uj) and /x(ti;) must be analytic in the upper half-plane 
due to causality requirements, eomg(w) and /Xemg(w) are 
not necessarily bound to causal behavior. In other words, 
it means that the analytical continuations of eemgC"^) and 
/^cmg(^) into the complex plane are not necessarily close 
to the continuations of e(cj) and /.t(w), and in this case 
they happen to be quite different. 

A possible way to proceed is to rely on the analytic 
properties of e(w) and write the Kramers-Kronig relation 



1 r^ 

e(w) = 1 + — P / dy 

IT^ Jo 



y~uj 



(68) 



where P stands for the Cauchy principal value, and con- 
sider also the analogous relation for fJ.{i-u). Taking the 
real part and evaluating it at an imaginary frequency i^, 
we obtain 



em = 1 + - / 

T^ Jo 



Imejy) 

'^yy'E2~, — 2' 



(69) 



and, using the fact that e{uj) ~ ecmg('^) [611, we have 



em ' 



IT 

2 

TT 



dyy 



Im Eo 



M 



e+y^ 



Im/icmg(2/) 

dyy g2 I 2 ■ 
C^ + r 



(70) 



FIG. 11: The plot of AP = P<^' - P'^' for different tem- 
peratures and models. Following our conventions throughout 
the paper, a positive force means attraction. The param- 
eters of P^^' are the same used in Fig. 8, in dimensional 
units they are Od = 1.32 x 10^'^rad/s, 70 = 5.48 x 10"rad/s 
(70 = for the plasma curve), fie = 4.7 x lO^^rad/s, ftm = 
8.7 x 10"rad/s, uj^ = 2.7 x lO^^rad/s, cjm = 2 x lO^^rad/s, 
7e/ = 5.5 X 10"rad/s, 7™ = 2.7 x 10"rad/s, f^a.i = 2.52 x 
lO^'^rad/s, LJ2.1 = 2.48 x 10"rad/s, 02,2 = 6.4 x lO^^rad/s, 
1.^2,2 = 3.8 X lO^^rad/s, ^2,3 = 1^2.3' = 1-9 x 10"rad/s, 
72,1 ~ 72,2 = 72,3 = 0, / = 0.1, and the parameters of P'^' 
are exactly the same except for Qe = 0. The inset shows the 
same plot on a different scale, since in the larger one it is not 
possible to see AP for large distances. 



In Fig. 10 we plot e(i^) and iiji f) u sing approximate 
values for TlCl polaritonic spheres [56[ embedded in vac- 
uum. We sec that e(i^) is overwhelmingly dominant over 
/i(i^), which in fact is hardly different from unity. As 
the insets show, this is basically due to a single strong 
resonance, around to — 0.317, that appears in eemg('^) but 
not in /.tomg(w)- From these results we conclude that, de- 
spite the fact that some magnetic activity is created by 
an array of polaritonic spheres, the Casimir force in this 
case is dictated by the electric part alone and therefore 
no repulsion seems possible. 



V. DISCUSSION 

A striking conflrmation of the magnetic influence on 
the Casimir force would be a measurement of repulsion 
between a metallic plate and a magnetodielectric one. 
This seems unlikely in light of the examples presented 
here, but a measured reduction in the attractive force 
might nevertheless be traced back to the magnetic prop- 
erties of a metamaterial. 



Let P'"'^' be the Casimir pressure between a gold half- 
space and a given metamaterial. If the magnetic prop- 
erties of the MM are "turned off" , keeping all other 
parameters the same, the pressure will change to some 
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FIG. 12: The difFerence in frequency shifts caused by the 
presence or absence of magnetic activity in the metamaterial. 
Everything is assumed to be at zero temperature. The MM 
is the same used in Fig. 11, and the parameters for the Rb 
atom are m = 1.45 x 10"^^ Kg, ao = 4.74 x 10"^^ cm^ and 
u>o = 2.54 X lO'^^ rad/s, with an unperturbed trap frequency 
of a;^ = 27r X 229 Hz. 



typical bulk-bulk measurements, which is desirable from 
the point of view of an effective medium approximation. 
The zero temperature Casimir-Polder potential between 
a ground state atom and a material half-space is [64| . 



22^-3 



f^cp(^) = ^^ / d^eam / dkk—- 



TE.TE 



{l^.k)-[l + 



„TM,TM 



{li^k) 



(71) 



where z is the distance between the atom and the half- 
space, K'i is defined just below dH]), and r^'^''^^ and 
j,TM,TM g^j.^, gjyQ]^ |-,y ([Y]) a{i^) is the dynamic atomic 

polarizability, which we assume is described reasonably 
well by the single-resonance expression 



c^m 



tto 



^+e/u^v 



(72) 



new value P^'^\ In order to check whether the differ- 
ence AP = P'^' — p(^) should be observable, we plot its 
computed value in Fig. 11 for zero and room tempera- 
tures, using both Drude and plasma models for the metal. 
The sensitivity of current experiments lies around 1 mPa, 
from which we conclude that detection of magnetic effects 
in our setup is currently possible up to d ~ 0.4/im. While 
this suggests a considerably large window for measure- 
ment, given that many experiments probe the 150 — 350 
/im range quite accurately, several things must be dealt 
with. First and foremost, we see that the difference be- 
tween the Drude and plasma predictions are consider- 
ably large (as compared to the magnetic effect) above 
0.6/im. This means that in order to ascribe changes in 
the Casimir force ambiguously to magnetic effects one 
has to know how to model metallic materials properly. 
In addition, at close distances like d < QAixm, the ef- 
fective medium approximation probably no longer holds, 
since the very structures that produce magnetic activ- 
ity (the metallic spheres in this example) are built on 
the scale of hundreds of nanometers or larger. These 
finite-size effects should bring significant corrections to 
the Casimir force, and must be considered in a more so- 
phisticated analysis. Finally, there are the imperfections 
of the materials themselves, like roughness, that at those 
distances play a non-negligible role. We conclude then 
that despite the fact that current experiments have in 
principle the sensitivity necessary to detect magnetic ef- 
fects, an actual measurement of such effects remains a 
challenging task. 



Casimir-Polder experiments [62l . |63[ also provide pos- 
sibilities for the detection of magnetic effects. These ex- 
periments are able to probe larger distances than the 



where a^ is the static polarizability and ujq is the dom- 
inant atomic transition. In one type of experiment [63[ 
the directly measured quantity is the frequency shift in 
the center of mass oscillation of a Bose-Einstein conden- 
sate: 



7(2) 



1 



2muj'. 



dlUcp{z')\.^^ 



(73) 



where m is the atomic mass and lOz is the unper- 
turbed (i.e., without Casimir-Polder forces) oscillation 
frequency. The reported sensitivity for 7 lies between 
10~^ and 10~^, setting the lower bound for the detec- 
tion of magnetic effects in the Casimir force. Let us then 
consider a Rb atom in front of the same MM used in 
the previous example, and compare the frequency shifts 
when its magnetic part is "turned on" and "off" . In Fig. 
12 we plot the difference A.'j(z) = "fnm{z) — jmiz), where 
7m (z) and jnmiz) are respectively the frequency shifts 
when the magnetic activity is present and absent. We 
see that in the best case scenario (sensitivity equal to 
10~^) the magnetic influence would be detectable up to 
around 2.5 /im; for larger distances the force is just too 
weak. 

As a final remark we note that, while Casimir repul- 
sion will likely be very difficult to observe with existing 
mctamaterials, the detection of magnetic effects through 
a slight reduction in the Casimir attraction is definitely 
possible. There are still some issues to be dealt with, like 
the assurance that magnetic activity is the main cause 
of force reduction, rather than some trivial effect like a 
reduced filling factor. With the consistent development 
of both Casimir measurements and MMs manufactured 
in the recent years, it is very reasonable to expect that a 
Casimir measurement of magnetic effects will be feasible 
in the near future. 
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